NONLINEAR STOCHASTIC BIASING OF GALAXIES AND DARK HALOS 
IN COSMOLOGICAL HYDRODYNAMIC SIMULATIONS 
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ABSTRACT 



We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and 
dark halos in spatially flat low-density CDM universe (Oq = 0.3, Aq = 0.7, /i = 0.7, 
and a% = 1) using cosmological hydrodynamic simulations. We identify galaxies by 
linking cold and dense gas particles which satisfy the Jeans criterion. We compare their 
biasing properties with the predictions of an analytic halo biasing model. Dark halos in 
our simulations exhibit reasonable agreement with the predictions only on scales larger 
than ~ 10/i^^Mpc, and on smaller scales the volume exclusion effect of halos due to 
their finite size becomes substantial. Interestingly the biasing properties of galaxies are 
better described by extrapolating the halo biasing model predictions. 
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The clustering amplitudes of galaxies are almost independent of the redshift between 
z = and 3 as reported in previous simulations. This in turn leads to a rapidly evolving 
biasing factor; we find that 6cov — 1 at redshift z ~ to 6cov — 3 — 4 at z = 3, where 
6cov is a biasing parameter defined from the linear regression of galaxy and dark matter 
density fields. Those values are consistent with the observed clustering of Lyman-break 
galaxies. 

We also find the clear dependence of galaxy biasing on their formation epoch; the 
distribution of old populations of galaxies tightly correlates with the underlying mass 
density field, while that of young populations is slightly more stochastic and anti-biased 
relative to dark matter. The amplitude of two-point correlation function of old popu- 
lations is about 3 times larger than that of the young populations. Furthermore, the 
old population of galaxies reside within massive dark halos while the young galaxies 
are preferentially formed in smaller dark halos. Assuming that the observed early and 
late-type galaxies correspond to the simulated old and young populations of galaxies, re- 
spectively, all of these segregations of galaxies arc consistent with observational ones for 
the early and late-type of galaxies such as the morphology-density relation of galaxies. 

Subject headings: galaxies: clustering - galaxies: formation - galaxies: halos - dark 
matter - cosmology: large-scale structure of universe - methods: numerical 

1. INTRODUCTION 

Ongoing galaxy redshift surveys, such as the Sloan Digital Sky Survey (SDSS) and the Two- 
Dcgrcc Field (2dF), aim at revealing the the large-scale structure of the universe with unprecedented 
precision. The gravitational instability is the main key process of the dark matter clustering, and 
this is now well understood from numerical simulations and several empirical theoretical models 
(Davis et al. 1985; Hamilton et al. 1991; Suto 1993; Mo & White 1996; Navarro, Frenk & White 
1997). In fact once the underlying cosmological models are specified, the two-point correlation 
functions of dark matter, which are the most conventional and widely used statistics describing the 
large-scale structure, can be fairly accurately predicted even with the redshift distortion and light- 
cone effects (Peacock & Dodds 1996; Suto et al. 1999; Suto, Magira, Sz Yamamoto 2000; Hamana, 
Colombi, k Suto 2001). 

On the other hand, it is widely believed that the distribution of galaxies is somewhat biased 
with respect to the underlying dark matter. For instance, Lyman-break galaxies at redshift z ~ 3 
(Steidel et al. 1998; Adelberger et al. 1998; Giavalisco et al. 1998) exhibit strong clustering and 
the galaxy biasing with respect to dark matter is time-dependent. Also the galaxy clustering is 
dependent on the galaxy morphology and environment (Dressier 1980; Postman & Cellar 1984; 
Loveday et al. 1995; Hermit et al. 1996; Dressier et al. 1997; Tegmark &: Bromley 1999) indicating 
the galaxy biasing is sensitive to many physical processes and thus stochastic. Clearly the rela- 
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tion between galaxy and dark matter clustering is far from simple and not yet fully understood 
either observationally or theoretically. This is the primary difficulty in properly interpreting the 
observational data of the upcoming large-scale redshift surveys. 

So far several models of galaxy biasing have been proposed adopting simplifying assumptions; 
Fry (1996) and Tcgmark &: Peebles (1998) discuss the evolution of biasing assuming that the 
number of galaxies does not change. Mo &: White (1996) present a model for the nonlinear biasing 
of virialized dark halos using the extended Press-Schechter formalism (Bond et al. 1991). Jing 

(1998) tested and improved the formula for the biasing of halo correlation functions originally 
proposed by Mo k. White (1996) using high-resolution N-hodj simulations. Dekel k. Lahav (1999) 
develop a fundamental framework to quantify the nonlinearity and stochasticity in galaxy biasing. 
Their formulation was subsequently applied to several numerical simulations (Blanton et al. 1999, 
2000; Somcrville et al. 2001). The biasing of dark halos is also investigated by Kravtsov & Klypin 

(1999) using high resolution N-hodj simulations. Note that their definition of dark halos is different 
from the conventional one used in the Press-Schechter formalism but rather close to dark matter 
cores (DM cores) in our analysis below. Recently, Taruya & Suto (2000; TS hereafter) proposed a 
first physical and analytical model for nonlinear and stochastic halo biasing combining the biasing 
model of Mo & White (1996) and the formation epoch distribution (Kitayama k. Suto 1996). 

More realistic approaches to galaxy biasing employ the state-of-the-art numerical simulations 
including the mcsh-bascd hydro dynamical simulations (Blanton et al. 1999, 2000; Cen & Ostrikcr 
2000) , and A^-body simulations combined with semi-analytic modeling of galaxy formation (Benson 
et al. 2000; Somerville et al. 2001). In what follows, we use the cosmological smoothed particle 
hydrodynamic (SPH) simulations (Yoshikawa, Jing k Suto 2000) of cold dark matter (CDM) uni- 
verse to examine the galaxy biasing. In particular, we focus on the comparison of the biasing 
characteristics of simulated objects (galaxies and dark halos) with the halo biasing model of TS. 
In addition, we investigate dependence of galaxy biasing properties on their formation history as 
an origin of galaxy morphology. Our simulation directly follows hydrodynamical and radiative pro- 
cesses to simulate galaxy formation, while the evolution of galaxies is not so properly modeled as 
those combined with a semi-analytic method of galaxy formation (Somerville et al. 2001). Due to 
the Lagrangian nature of the SPH technique, the spatial resolution of our simulations is better than 
those in the mesh-hydro simulations by Blanton et al. (1999, 2000) and we can resolve galaxies as 
distinct and isolated objects, while their treatment of the thermal process and the metal enrichment 
is more realistic. Thus our method is complementary to those previous investigations with different 
approaches. 

The rest of the paper is organized as follows. In §2, we describe the detail of our numerical 
simulation and the procedures to identify galaxies, dark matter cores (DM cores), and dark halos. 
In §3 we present a brief summary of the biasing description following TS, and compare several 
properties of the biasing in the one-point statistics of galaxies and dark halos. Then we discuss 
the biasing in terms of their two-point correlation functions. Section 4 examines the dependence 
of galaxy biasing on their formation history. Finally, we summarize our major findings in §5. 
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2. METHODOLOGY 

2.1. Cosmological SPH simulation 

Our numerical simulation code is a hybrid of Particle-Particle Particle-Mesh (P^M) A^-body 
Poisson solver (Hockney & Eastwood 1981) and smoothed particle hydrodynamics (SPH) algorithm 
(Yoshikawa, Jing &: Suto 2000). The simulation presented in this paper adopts A^dm = 128^ dark 
matter particles and the same number of gas particles for SPH. We use the spline (S2) functional 
form for gravitational softening (Hockney & Eastwood 1981) and the softening length is set to 

1/3 

Egrav = -^box/(10A/Q]y[) and kept constant in comoving coordinates, where -Lbox is the comoving 
size of the simulation box. We set the minimum of SPH smoothing length to /imin = egrav/4 and 
adopt the ideal gas equation of state with an adiabatic index 7 = 5/3. The effect of radiative 
cooling is included adopting the metallicity of [Fe/H]= —0.5. We use the cooling rate described 
in Sutherland & Dopita (1993). Thacker et al. (2000) reported that artificial over-cooling occurs 
under the presence of radiative cooling in SPH simulations due to overestimate of hot gas density in 
the vicinity of cooled gas clumps due to the smoothing scheme of SPH algorithm. In order to avoid 
this numerical artifact, we implement a modification of SPH algorithm, "cold gas decoupling", 
following Pearce et al. (1999). The detail of this prescription is presented in the next subsection. 

We consider a spatially-fiat low-density CDM (LCDM) universe with JIq = 0.3, Aq = 0.7, ag = 
1.0 and h = 0.7, where is the mean mass density parameter, Aq the dimensionless cosmological 
constant, cxg the rms density fluctuation on a scale of 8hr^ Mpc and h the Hubble constant in 
units of 100 km-s^^-Mpc^^. This particular model satisfies both the COi?i? normalization (Bunn &: 
White 1997) and the abundance of X-ray clusters of galaxies (Kitayama & Suto 1997). We assume 
the mean baryon mass density parameter to be Qb = 0.015h~^ (Copi et al. 1995). The simulation 
is carried out in a periodic cube of (75/i~^Mpc)^, with the gas and dark matter mass per particle 
being 2.4 x 10^ M0 and 2.2 x 10^^ Mq, respectively. The initial condition is created at 2; = 25 using 
the COSMICS package (Bertschinger 1995), which is evolved up to 2: = 0. 



2.2. Cold gas decoupling and identification of galaxies 

In order to avoid the numerical over-cooling of gas particles mentioned above, we decouple 
cold gas particles which satisfy the following Jeans condition (Yoshikawa, Jing & Suto 2000): 

hsPH > I ' (1) 

where /isph is the smoothing length of gas particles, Cg the sound speed, G the gravitational constant 
and pgas the gas density of gas particles. Except for the fact that these cold gas particles are ignored 
in computing the gas density of hot gas particles, all the other SPH interactions are left unchanged. 
This decoupling scheme is a phenomenological treatment of multi-phase gas dynamics, and should 
be interpreted as an approximate prescription of galaxy formation. 
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Galaxies in our simulations are identified as clumps of cold and dense gas particles which 
satisfy the criterion (1) and 

Pgas > 10^Pb(2), (2) 

where p\,{z) is the mean baryon density at redshift z. Figure 1 shows the scatter plots of gas 
particles in density - temperature plane. The blue points indicate the cold and dense gas particles 
satisfying the criteria (1) and (2), the diffuse cold gas particles which satisfy (1) and pgas < 10^ Ph{z) 
are plotted in green, and the other hot gas particles are represented in red. This indicates that 
the above criteria for the galaxy particles properly segregate the cold and dense gas particles. We 
group these particles using fricnd-of-friend (FOF) algorithm (Davis et al. 1985) with linking length 
bg = 0.0164(1 + z)l and identify the resulting clumps as "galaxies", where I = Lbox/^DM 
comoving mean particle separation. The proper choice of the linking length is not clear and we 
simply adopt the value of Pearce et al. (1999) here. In this paper, we only consider galaxies with 
mass greater than Mg^min = 10^^ Mq, which is equivalently 40 times of each gas particle mass and 
close to a nominal mass resolution of baryonic matter^. As noted in § 2.3, the mass functions 
of simulated galaxies are roughly consistent with those from semi-analytic modeling of galaxy 
formation, which justifies our galaxy criteria empirically to some extent. We show the number of 
galaxies identified in our simulation and the adopted linking length in Table 1. 



2.3. Identification of dark halos and dark matter cores 

The FOF algorithm is also applied in identifying dark halos. The linking length 6h for dark 
halos is set to satisfy the equation 



187r2 V0.2r' ' 

where Ac{z) is the mean over-density of spherically virialized objects formed at redshift z. We 
compute Ac{z) at each redshift using a fitting formula by Kitayama & Suto (1996). At z = 0, for 
instance, Ac = 335 and bh = 0.164 T. 

We also identify the surviving high-density substructures in dark halos, which we call DM 
cores. Identification of substructures in dark halos is a technically challenging problem and several 
objective methods have been proposed so far (Gelb & Bertschinger 1994; Eisenstein & Hut 1998; 
Klypin et al. 1999; Springel et al. 2000). In order to identify DM cores in our simulation, we adopt 
the hierarchical FOF (HFOF) method (Gottlober et al. 1999). In HFOF method, we apply the 
conventional FOF method with a set of different linking length be- be = ^max/4, Zmax/2, and Imax, 
where Zmax is the maximum linking length. For each linking length, gravitationally bound groups 
with more than 20 particles are identified as DM cores. The maximum linking length is set to 

Zmax = 0.05r. 



^SPH gas density is smoothed over about 30 nearest neighbor gas particles. 
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In this paper, we consider the dark halos with their mass greater than 10^'^ Mq (~ f^o/^b x 
-^g,min) and DM cores with more than 20 dark matter particles (equivalently 4.3 x W^^Mq). In 
Table 1, the number of identified objects and adopted linking length are also shown. Figures 2 and 
3 show the distribution of dark matter particles, gas particles, dark halos and galaxies at z = 
and z = 2. At z = galaxies are more strongly clustered than dark halos, while at z = 2 those two 
objects show similar spatial distribution. 

Figures 4 and 5 show close-up snapshots of the most massive cluster at ^ = with mass 
M ~ 8 X 10^^ Mq and a relatively poor cluster with M ~ 10^^ Mq, respectively. In each figure, 
upper panels depict the distribution of dark matter and gas particles, and the distributions of DM 
cores and dense cold gas particles which satisfy the condition (1) and (2) are shown in lower panels. 
Circles in lower panels indicate the positions of galaxies identified in our simulation. We can see 
that for the richer cluster, the distribution of DM cores is relatively in good agreement with that 
of galaxies except for the cluster center, where the tidal radius is much shorter than our numerical 
resolution. On the other hand, galaxies or cold gas clumps in the smaller cluster are not necessarily 
hosted by DM cores. This is probably because DM cores in our simulation significantly suff^er from 
the artificial overmerging, which is severer for poorer dark halos due to small number of particles, 
while galaxies represented by dissipative gas particles are less affected by this overmerging. This 
is why DM cores at higher redshift are much less abundant than galaxies and dark halos (see 
Table 1). This problem is intrinsically related to the question of whether substructures within 
dark halos identified in high-resolution A''-body simulations (Klypin et al. 1999; Colm et al. 1999) 
really correspond to the real galaxies. Unfortunately the resolution of our current simulations is 
not sufficiently good to answer this issue in a reliably manner, but we hope to revisit this with 
another SPH run with N = 256^ particles (Yoshikawa, Jing Sz Suto, in preparation). 

Figure 6 shows mass function of dark halos and galaxies at z = and 2. We find that the mass 
function of simulated dark halos (upper panels) agrees better to the fitting formula of Jenkins et al. 
(2001) {dashed lines) than that of Press & Schechter (1974) (solid lines). Galaxy mass functions in 
our simulations (lower panels) are roughly consistent with those from other SPH simulations and 
semi-analytic models (Benson et al. 2001), but slightly less abundant at Mgaiaxy ^ lO^-'^M© due to 
limited mass resolution. 

3. BIASING PROPERTIES OF GALAXIES AND DARK HALOS 

The most natural form of galaxy biasing is the relation between over-density fields of galaxies 
Sg and dark matter S^i- In this section, we compute the density fields of galaxies and dark halos 
from our simulation, and study their statistical properties and evolution. 
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3.1. Formulation and computation of biasing parameters 

A biasing scheme relates the density field of dark matter with those of galaxies and dark halos, 
which are defined for a given smoothing scale Rg as 

C>m(X,-Ks) = z 1, (4) 

P 

<5g(x,i?,) = !^i^^-l, (5) 

(5h(x, Rs) = — ^ - 1, (6) 

nh 

where p{x,Rs), ng(x, and ni^{x,Rs) denote the mass density, and galaxy and halo number 
densities smoothed over the top-hat window radius Rg, and the over-bar (~) indicates the mean over 
the entire universe. We briefly summarize several parameters describing the nonlinear stochastic 
nature of biasing introduced by TS for later convenience. 

The joint probability distribution function (PDF), P{6m,Si), characterizes the statistical prop- 
erties of (5in and Si, where the subscript i indicates two different objects; "g" for galaxies and "h" 
for dark halos. By definition, Sjn and Si have zero mean and their variances are related to the joint 
PDF as 

f^m = (<^m) = J J P{S^,Si)SldSradSi (7) 

and 

= (<^l) = j j Pi^m, Si)Sf dS^ dSi, (8) 

where the bracket (■ ■ ■) denotes the joint average over Si and S-ai- The statistical relation between 
Si and Sra is described by the conditional PDF, P{Si\S-ai)- The conditional mean of Si, Si{S-ai), for a 
given S-ai is then calculated from 



Si{Sm) = j S^P{S^\S,,)dSi , (9) 
yielding the following biasing parameter (TS): 

, _ {Si{Sra)Sra) _ jSiS^) 

"cov.j — 9 — 9 ■ \^^) 

The nonlinearity of the biasing is quantified by 

2 _ (0(^7) <m , 

{SiSr.)^ {SiS^)^ ' ^ ' 

which vanishes only when the biasing is linear (i.e., the ratio Si/Sjn is independent of S^) and is 
positive otherwise. Similarly the stochasticity of the biasing is characterized by 

2 _ {si)m-m _ <wf-m ..9^ 
(^,,5^)2 {SiS^}-^ • ^ ' 
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This parameter vanishes for the deterministic bias where 5i = 5i{5-ai)- In terms of the above biasing 
parameters, a somewhat more conventional biasing coefficient 6var,i = (^i/(^ra is written as 

^var,i = bcov,i (1 + ^n\,i + ^scatt.i) ^ • (l^) 

Finally the correlation coefficient r-corr.i (Dekel k. Lahav 1999) is given by 

_ {5i5^) 1 



(14) 



We compute the biasing parameters &cov,i, ^var,i, eni,i) Cscatt.i and rcorr,* each for dark halos and 
galaxies with smoothing scales R = 4/i^-'^Mpc, B/i^-'^Mpc and 12/?,^^Mpc. Wc obtain many pairs 
of the values ((5i(x, i?s), (5in(x, i?s)) for randomly selected points x in the simulation volume and 
evaluate the biasing parameters using equations (7) - (14) by replacing the joint averages (• • •) 
with averages over all selected points. The number of randomly selected points is 1000 for the 
top-hat smoothing scale Rg = 12^-^Mpc, 5000 for Rg = 8/i~^Mpc and 30000 for Rg = 4^-^Mpc. 
Since our simulation volume is 75/i~^Mpc per side, most of the selected sampling points are not 
fully independent. Nevertheless we decided to make over-sampling in evaluating the mean and the 
variance of the density fields. Thus our quoted error-bars below may rather correspond to those in 
a bootstrap resampling method. 



3.2. Comparison of biasing of galaxies and dark halos 

Figure 7 shows the joint distribution of and 6g with mass density field dm at redshift z = 0, 
1 and 2 smoothed over Rg = 12^~^Mpc {Upper panels) and 4/i~^Mpc {Lower panels). We plot the 
conditional mean relation Si{6in) from our simulation results {solid lines) and from the theoretical 
prediction of halo biasing by TS {dashed lines). In computing theoretical predictions, we adjust 
the range of dark halo mass as our simulated dark halos (Table 1). 

Consider first the results for dark halos. For a given smoothing scale, the simulated halos 
exhibit positive biasing for relatively small in agreement with the predictions. On the other 
hand, they tend to be underpopulated for large Sm, or anti-biased. This is mainly due to the 
exclusion effect of dark halos due to their finite volume size as previously discussed in Taruya et 
al. (2001) using purely N-body simulations. The theoretical model of TS does not take account of 
this effect, and thus the discrepancy between the predictions and the simulations becomes more 
substantial for smaller Rg and/or at lower z as expected. 

Since our identified galaxies have smaller spatial extent than the halos, the exclusion effect is 

not so serious. This is clearly illustrated in lower panels in Figure 7. In fact they seem to show 
much better agreement with the TS predictions despite the fact that the models are formally valid 
only for dark halos defined according to the Press-Schechter manner. 
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A more careful look at the results for galaxies, however, reveals that 5g(5m)/5m decreases 
slightly at larg er 6ixi especially for smaller smoothing scale — Ah "'^Mpc. While this tendency 
may be partially explained by their volume exclusion effect, their typical sizes seem to be sufficiently 
small to account for this. Rather, we consider two possible origins of this tendency. One is the 
suppression of galaxy formation at very high temperature and thus high density regions, as pointed 
out in Blanton et al. (1999, 2000). Figure 8 shows the dependence of galaxy overdensity on the 
surrounding gas temperature separately for galaxies with different formation redshifts (see §4 for 
details), and supports this interpretation; the ratio, (1 + Sg)/{1 + 6m), is anti-correlated with the 
surrounding gas temperature. Comparing the left and right panels in Figure 8 indicates that 
the anti-correlation with gas temperature is much stronger for galaxies which form relatively late 
{z{ < 1.7). On the other hand, those formed earlier show very weak, at most, anti-correlation, 
which is natural because they should have collapsed and formed much before the surrounding gas 
acquires the current high temperature. Another possibility is that there is an intrinsic difference in 
formation epoch of galaxies between over- and under-dense regions. Since, in hierarchical formation 
scenarios, objects in over-dense regions tend to form earlier than those in under-dense regions, it 
is expected that young galaxies with Z{ < 1.7 form relatively lower-dense thus low-temperature 
regions, which is also consistent with Figure 8. The similar analysis for DM cores will distinguish 
these two possibilities. Although we notice that the same correlation exists even for DM cores, we 
suspect that this is mainly due to the artificial overmerging effect as we discussed in §2.3., and 
will revisit this topic with another simulation with higher-resolution (Yoshikawa, Jing &; Suto, in 
preparation). 

Incidentally, in order to check the dependence of the simulated galaxy biasing on the lower 
mass limit of our criteria, Mgaiaxy > lO^-'^AfQ (or equivalently iVgas > 40), we construct another 
set of galaxy sample adopting higher mass cutoff A^gas > 80, and compare their biasing properties. 
We find that the joint probability distribution of and Sg for the galaxy sample selected with 
A^gas > 80 does not significantly change from those of the original galaxy sample. 



3.3. Stochasticity and nonlinearity in biasing of galaxies and dark halos 

The stochasticity and nonlinearity in galaxy and halo biasing are clearly identified in Figures 7. 
For more quantitative discussion, we plot in Figure 9 the evolution of their biasing parameters 6cov) 
rcorr, Escatt and ecorr for three different smoothing radii. 

Consider first 6cov This biasing parameter exhibits strong time-dependence; the biasing is 
stronger in the past. This is consistent with analytic biasing models (Mo &; White 1996; Taruya 
& Suto 2000), previous numerical simulations (Kravtsov &: Klypin 1999; Somerville et al. 2001; 
Pearce et al. 1999) and in fact explains the recent observations of Lyman break galaxies (Giavalisco 
et al. 1998; Adelberger et al. 1998) using the halo biasing model (Mo & White 1996; Jing & Suto 
1998). On the other hand, the scale-dependence of bcov is very weak as in the case of the biasing 
parameter defined through the two-point correlation function (see the next subsection). 
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Both escatt and Cni in our simulated catalogues are somewhat smaller than the TS prediction, 
but their qualitative behavior is consistent with the model; larger on small scales and almost 
independent of z. The biasing becomes linear and deterministic relation and also the volume 
exclusion is less effective for larger smoothing scales. The current degree of the stochasticity and 
nonlinearity does not hardly affect the amplitude of clustering (see eq.[13]), but the topology of the 
isodensity contours is sensitive to the nonlinearity even at this level (Hikage, Taruya h Suto 2001). 

It is interesting to notice that the biasing parameters for galaxies show similar behavior and 
are closer to the predicted behavior. In addition, all biasing parameters for dark halos and galaxies 
behave very similarly at high redshifts ^; c± 2 — 3. This indicates that the spatial distribution of 
galaxies and dark halos are statistically similar, and can be understood by the fact that we have 
one-to-one correspondence between dark halos and galaxies at z ~ 2 — 3 as shown below. 

Figure 9 also shows that the evolution of biasing is almost independent on the lower mass limit 
of the galaxies. This might be interpreted as our simulated galaxy sample is nearly complete for 
the present purpose. 

Figure 10 shows the number of member galaxies which reside within the virial radius of their 
hosting dark halos {upper panels) and the distribution of their mass ratios (lower panels) at redshift 
z = 0, 2 and 3. Solid and dashed lines in lower panels indicate the cosmic mean baryon fraction 
Ob/^o and resolution limit of galaxy mass (Afgaiaxy = ^0^^ Mq), respectively. One can see that 
most dark halos at 2; = 3 host only one galaxy, explicitly justifying the empirical assumption of 
one-to-one correspondence between dark halos and Lyman-break galaxies around z = 3 in previous 
studies (Jing & Suto 1998; Steidel et al. 1998). The subsequent evolution of dark halos involves 
several merger processes, and thus dark halos at lower redshifts tend to host multiple member 
galaxies. 



The previous subsections discuss only the biasing parameters defined from the one-point statis- 
tics. In this subsection, we turn to a more conventional biasing parameter defined through the 
two-point statistics: 



where ^ii{r) and ^nim{f) are two-point correlation functions of objects i and of dark matter, re- 
spectively. While the above biasing parameter is ill-defined where either ^^(r) or ^mm(?") becomes 
negative, it is not the case at clustering scales of interest (< lOft-^^Mpc). The relation of one-point 
and two-point biasing parameters is also investigated in detail by Taruya et al. (2001) for density 
peaks and dark halos. 



3.4. Biasing in terms of the two-point correlation function 




(15) 



Figure 11 shows two-point correlation functions of dark matter, galaxies, dark halos and DM 
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cores {upper and middle panels), and the profiles of biasing parameters b^{r) for those objects 
(lower panels) at z = 0, 1 and 2. In the upper panels, wc show the correlation functions of DM 
cores identified with two different maximum linking length; l^ax = 0.05 as presented in § 2.3 and 
^max = bh/2. Correlation functions of DM cores identified with /max = 0.05 are similar to those of 
galaxies. On the other hand, those identified with /max = &h/2 exhibit much weaker correlation, 
and are rather similar to those of dark halos. This is due to the fact that HFOF algorithm with 
larger Zmax tends to pick up lower mass halos which are poorly resolved in our numerical resolution. 

The correlation functions of galaxies are almost unchanged with redshift, and that of dark 
halos only slightly evolves between z = and 2. By contrast, the amplitude of the dark matter 
correlation function evolves rapidly by factor of ~ 10 from z = 2 to z = 0. The biasing parameter 
b^^g is larger at a higher redshift, for example, b^^g ~ 2-2.5 at z = 2. These results are consistent 
with the numerical studies by Bagla (1998), Colin et al. (1999) and Pearce et al. (1999) and also 
qualitatively explains the clustering of Lyman-break galaxies (Giavalisco et al. 1998). The biasing 
parameter b^^h for dark halos is systematically lower than that of galaxies and DM cores again due 
to the volume exclusion effect. At z = 0, galaxies and DM cores are slightly anti-biased relative 
to dark matter at r ~ l/i~^Mpc, which is also consistent with previous numerical simulations 
(Pearce et al. 1999; Colin et al. 1999; Benson et al. 2000; Somerville et al. 2001) and also with the 
observational results from the Las Campanas Redshift Survey (Jing, Mo & Borner 1998). In lower 
panels, we also plot the one-point biasing parameter 6var,i ^ Ci/fTm at r = Rg for comparison. In 
general we find that b^^i is very close to 6var,« at z ~ 0, but systematically lower than 6var,i at higher 
redshifts. 



4. THE FORMATION EPOCH AS AN ORIGIN OF THE MORPHOLOGICAL 

TYPE OF GALAXIES 

It is fairly established that there exists a certain correlation between the morphology of galaxies 
and their star formation history; early-type galaxies form via initial star bursts at high redshifts 
while late-type galaxies experience continuous and relatively mild star formation history (Roberts & 
Haynes 1994; Kennicutt 1998). This implies that the galaxy morphology is empirically related to its 
formation epoch. On the basis of this interpretation, one can examine the morphology-dependent 
clustering of galaxies by classifying our simulated galaxies according to their formation epoch. 

We have fifty outputs of all simulation particles at different redshifts between z = 9 and 0. For 
each galaxy identified at z = 0. wc define its formation redshift Zf by the epoch when half of its 
cooled gas particles satisfy the criteria (1) and (2). Roughly speaking, Zf corresponds to the median 
formation redshift of stars in the present-day galaxies. We divide all simulated galaxies at z = 
into two populations (the young population with Z{ < 1.7 and the old population with Zf > 1.7) so 
as to approximate the observed number ratio of 3/1 for late-type and early- type galaxies (Loveday 
et al. 1995). 
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Figure 12 shows the joint probabihty distribution of Sm. and Sg respectively for the old {left 
panel) and young {right panel) populations. They exhibit clear difference in their clustering prop- 
erties. Their biasing parameters are Ug = 1.73(1.06), 6var,g = 1.51(0.93) and rcorr,g = 

0.95(0.88) 

for the old (young) population. These results qualitatively agree with Blanton et al. (1999), and 
Somerville et al. (2001) also showed a similar result that red galaxies are biased compared to the 
overall population and blue ones are anti-biased, where galaxies with color B — V > 0.8 are defined 
as red galaxies and the remainder as bhic ones. The dashed lines in both panels indicate the TS pre- 
dictions of the mean biasing for dark halos restricting the formation epoch as Zf > 1.7 and Zf < 1.7, 
respectively; the old population shows excellent agreement with the halo biasing prediction while 
the young population behaves rather differently. This indicates that early-type galaxies preferen- 
tially reside in the center of the massive halos almost in a one-to-one manner and that late-type 
galaxies avoid the dense environment, which is consistent with the observed morphology-density 
relation (Dressier 1980; Postman & Cellar 1984; Dressier et al. 1997). 

This interpretation is directly confirmed in Figure 13. Massive halos have a larger fraction 
of the old population of galaxies, while the young population of galaxies mainly reside in smaller 
halos. This segregation may be understood by the same mechanisms of anti-bias of galaxies at 
high density regions. As discussed in § 3.2, due to the suppression of galaxy formation in high 
temperature regions at lower redshift, and/or a different formation epoch for over- and under-dense 
regions, we have the deficiency of the young population of galaxies within massive dark halos at 
z = 0, and galaxies formed at high redshift trapped within the gravitational potential of dark halos 
gradually tend to trace the distribution of underlying dark matter. 

The difference of the clustering amplitude can be also quantified by their two-point correlation 
functions at 2; = as plotted in Figure 14. The old population indeed clusters more strongly than 
the mass, and the young population is anti-biased. The relative bias between the two populations 
ft^'^g = A/Cold/'^young ranges 1.5 and 2 for 1/i^^Mpc < r < 20h^^Mpc, where ^young and ^old are 
the two-point correlation functions of the young and old populations. Again this is completely 
consistent with the observational indications that the clustering of early- type galaxies is stronger 
than that of late-type galaxies by a factor of 3-4 in terms of the amplitude of two-point correlation 
functions (Loveday et al. 1995; Hermit et al. 1996). 

All the above results suggest that the old and young populations of galaxies in our simulations 
may be interpreted as the early-type and late-type galaxies in the present universe, and that the 
formation epoch and the hydro dynamical environment play the important role in determining the 
morphology of galaxies. We note here that the above result is fully consistent with the recent 
analysis of the IRAS PSCz galaxy sample by Jing, Borner, & Suto (2001), who found a strong 
anti-bias of the IRAS-selected galaxies (and thus mainly late- types). The degree of the detected 
bias is accounted for by the phenomenological cluster-underweight bias model (Jing, Mo & Borner 
1998), and also by the semi-analytic modeling of galaxy formation which assumes that the galaxy 
morphology is determined by the frequency of the major merger of halos (Kauffmann, Chalrot & 
White 1996; Kauffmann, Nusser & Steinmetz 1997; Kauffmann et al. 1999). 
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5. CONCLUSIONS AND DISCUSSION 

Using a cosmological SPH simulation, we directly simulate the formation of galaxies via radia- 
tive cooling of baryonic component and identify galaxies as isolated and distinct groups of cold gas 
particles. We calculated the biasing of galaxies and dark halos, and in particular, compared their 
properties with the theoretical prediction of the halo biasing model proposed by TS. 

Our major findings are summarized as follows; 

(1) The clustering of dark halos suffers from the the volume exclusion effect due to their finite 
size, especially at small scales. On the other hand, the halo biasing model by TS can reasonably 
account for the clustering of "galaxies" at large scales. At smaller scales, however, galaxies are 
anti-biased relative to dark matter at high density and thus high temperature environment. 

(2) The biasing parameters are strongly time-dependent. At z ~ 3, our galaxies exhibit strong 
biasing; 5cov,g — 3-4 and b^^g ~ 3, which is consistent with the observed clustering of Lyman-break 
galaxies (Adelberger et al. 1998; Giavalisco et al. 1998; Steidel et al. 1998). 

(3) The formation epoch Zf is the major parameter in determining the morphological type 
of galaxies. In our specific example, galaxies identified at z = with Zf > 1.7 and Z{ < 1.7 
can be roughly regarded as early-type and late-type galaxies, respectively. The former tightly 
correlates with the massive host halos and shows stronger clustering, while the latter is anti-biased 
and more stochastic. These suggest that biasing properties of galaxies, identified by different 
photometric bands or color selections, should be significantly different, which should be kept in 
mind in comparing the galaxy clustering from different galaxy catalogues. 

Our current definition of galaxies in simulation data is admittedly rather phenomenological. 
Apparently more observationally oriented classification of galaxies, for example using color or mag- 
nitude of galaxies, is necessary for direct comparison with observations. We plan to implement 
more realistic prescriptions of galaxy formation and evolution including star formation, feedback 
and UV background heating in due course. Nevertheless it is quite encouraging that even a simple 
scheme described here explains the major properties of galaxy clustering in the universe. 
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Table 1. Number and mass range of identified objects and adopted linking length for FOF 

algorithm. 



redshift 


dark halo 


galaxy 


DM core 


0.0 


1797 (0.164 [) 


1604 (0.01640 


1525 (0.05 




IO^^Mq ~ 8.6 X 10^^ M© 


lO^^Mo ~ 9.5 X IO^^Mq 


4.3 X lO^^Mo - 2.0 X 10^^ Mo 


0.5 


2105 (0.1840 


1936 (0.0246 


1721 (0.05 




10^2 ~ 3.5 X 10^4 Mo 


lO^^Mo ~ 6.8 X 10^2 j^^o 


4.3 X lO^^Mo ~ 6.3 X W^Mq 


1.0 


2201 (0.192 


1861 (0.0328 


1543 (0.05 




IO^^Mq ~ 2.3 X 10^^ Mo 


IO^^Mq - 3.6 X 


4.3 X 10" Mq ~ 2.9 X W^^Mq 


2.0 


1859 (0.197[) 


1360 (0.0492 


765 (0.05 




10^2 ~ 8.2 X 10^^ Mo 


lO^^Mo ~ 2.0 X lO^^Mo 


4.3 X lO^Mo ~ 9.8 X IO^^j^^q 


3.0 


1165 (0.199 


996 (0.0656 


278 (0.05 




10^2 ~ 3.8 X 10^3 Mo 


10^^ Mo ~ 1.4 X 10^2 Mo 


4.3 X lO^Mo ~ 4.2 X IO^^j^^q 
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Fig. 1. — Scatter plot of gas particles at z = in logpgas^ogT plane. The cold and dense gas 
particles which satisfy the criterion (1) and pgas > 10^ Phi^) are indicated by blue points, the diffuse 
cold gas with pgas > 10^ phi^) are by green points and the others are by red points. 
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Fig. 2. — Distribution of gas particles, dark matter particles, galaxies and dark halos in the volume 
of 75h~^ X 75h~^ x 30^~^Mpc^ model at 2; = 0. Upper-rightigas particles; Upper-left: dark matter 
particles; Lower-right: galaxies; Lower-left: DM cores 



Fig. 3. — Same as Figure 2 but for z = 2. 



Fig. 4. — Snapshots of the most massive cluster (M ~ 8 x 10^"^ M©) in the simulation at z = 0. 
Upperleft: dark matter; Upper-right: gas; Lower-left: DM cores; Lower-right: cold gas. Circles in 
lower panels indicate the positions of galaxies identified according to our criteria. The comoving 
size of the box is 6.25/i~^Mpc per side. 



Fig. 5. — Same as Figure 4 but for a poorer cluster with M ~ 10 Mq. 
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Fig. 6. — Mass functions of dark halos and galaxies at redshift z = and z = 2. Solid lines 
in upper panels are theoretical predictions of Press-Schechter mass function and dashed hnes are 
fitting formula by Jenkins et al. (2001). 
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Fig. 7. — Joint probability distributions of overdensity fields for dark halos and galaxies with dark 
matter overdensity smoothed over Rg = 12/i~^Mpc (Upper panels) and Rg = Ah^^Mpc (Lower 
panels) at redshift z = 0, 1 and 2. Solid lines indicate the conditional mean 6i(6m) for each object. 
Dashed lines in each panel depict the theoretical prediction of conditional mean by Taruya Sz Suto 
(2000). 
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Fig. 8. — Relation between environmental temperature and the ratios of dark matter overdensity 
with that of galaxies with Zf > 1.7 {left panel) and Zf < 1.7 {right panel) in the high density regions 
(Jm > 5). Each point corresponds to the randomly selected point in the simulation box. The 
smoothing scale is set to Rg = 4/1" ^Mpc. 
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Fig. 9. — Evolution of biasing parameters bcov, ?'corr, Cscatt and eni for galaxies (solid lines), dark 
halos (dotted lines), and the theoretical predictions by Taruya & Suto (2000) (long dashed lines). 
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Fig. 10. — Number of galaxies which reside in a dark halo {upper panels) and distributions of mass 
ratio between galaxies and dark halos which host them [lower panels) at rcdshift z = and 2. Solid 
and dashed lines in lower panels indicate the mean baryon fraction (fib/^^o) and resolution limit 
constrained by the minimum mass of galaxies in their definition (see §2.2), respectively. 
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Fig. 11. — Upper panels show two-point correlation functions for dark matter and DM cores at 
redshift z = 0, 1 and 2. Middle panels for those of dark matter, galaxies and dark halos. In 
lower panels, the profiles of biasing parameter 6^(r) for dark halos, DM cores and galaxies at 
z = 0.0, 1 and 2 are shown. In lower panels, we also plot the parameter 6var on the smoothing scale 
Rs = 4/i~^Mpc, 8/i~^Mpc and 12/i~^Mpc at r = Rg for each kind of objects by different symbols. 
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Fig. 12. — Joint probability distributions of density fields of dark matter and galaxies with different 
formation epochs on the scale of Rs = 8/i~^Mpc. Left panel is for galaxies with Zi > 1.7 and right 
panel for ones with Z{ < 1.7. Solid lines indicate the simulated mean relations. For comparison, 
the predictions of mean biasing for dark halos with their formation redshift greater and less than 
1.7. are shown in left and right panels, respectively. 
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Fig. 13. — Same as Figure 10 except for old (left) and young {right) populations of galaxies at 
z = 0. 
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Fig. 14. — Two-point correlation functions for the old and young populations of galaxies at z = 
as well as that of dark matter distribution. The profiles of bias parameters b^{r) for both of the 
two populations are also shown in the lower panel. 
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